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An  analytic  study  of  the  general  stability  of  thin 
circular  cylindrical  shells  stiffened  by  a  system  of 
inclined  stiff eners  has  been  carried  out.   The  cylinder  is 
internally  pressurized  and  is  under  three  types  of  load- 
ings, namely,  axial  compression,  bending,  and  transverse 
shear  load.   The  stiffeners  are  not  very  close  but  are 
discretely  located  and  the  eccentricity  of  stiffeners  is 
disregarded.   The  method  of  the  solution  is  carried  out 
by  the  use  of  nonlinear  large  deflection  theory  and  the 
effects  of  initial  imperfections  in  the  strain-displace- 
ment equations  are  considered.   The  Ritz-method  is  used 

xi 


to  find  the  governing  equations  of  instability  of  stiffened 
shells  siobject  to  the  three  types  of  loading. 

Numerical  examples  for  cylinders  with  inclined 
stiff eners  at  different  angles  are  worked  out.   The  minimxim 
strength  of  shell  versus  different  rigidities  of  stiffener 
are  also  calculated.   The  relation  between  the  minimum 
strength  of  pressurized  shell  and  internal  pressure  is 
plotted.   The  computation  was  carried  out  on  the  IBM  709. 
Some  tests  have  been  made  on  the  shell  with  inclined 
stiff eners  at  angles  of  /  =  30°,  45°,  and  60°  under  axial 
compression  and  bending  loads.   Some  of  the  experimental 
data  are  plotted  to  compare  with  the  analytical  results. 
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CHAPTER  I 


INTRODUCTION 


Interest  in  shell  buckling  phenomena  dates  back  to 
1858  when  Fairbairn  {_!)   performed  his  experiments  on  the 
buckling  of  cylinders  under  uniform  external  pressure.   In 
1888,  the  first  theoretical  treatment  of  the  problem  was 
published  by  Bryan  (_2)  and  Love  (_3)  ,  and  since  that  time, 
many  investigators  have  been  attracted  to  the  various  prob- 
lems of  thin  shell  instability.   The  equations  for  thin 
shells  were  discussed  by  Timoshenko  (4.)  in  his  well-known 
book.   The  study  of  stability  of  isotropic  cylindrical 
shells  was  advanced  by  Southwell  {5,)  .   In  1933,  a  set  of 
equilibrium  equations  for  cylindrical  shells  subject  to 
torsion  was  derived  by  Donnell  {6)    based  on  the  theory  of 
small  deflections.   A  simplified  method  of  elastic-stabili- 
ty analysis  for  thin  cylindrical  shells  was  also  discussed 
by  Batdrof  Q) . 


if 
Underlined  numbers  in  parentheses  refer  to  entries 

in  the  references. 


Then  Suer  and  Harris  (_8)  applied  Donnell's  eighth- 
order  linear  equilibrium  equation  to  the  problem  of  the 
stability  of  cylinders  under  combined  torsion  and  hydro- 
static pressure.   They  obtained  a  solution  with  Galerkin's 
method  based  on  an  assumed  radial  deflection  function  in 
the  form  of  an  infinite  trigonometric  series. 

In  1961,  Seide  (j9)  had  presented  an  analysis  of 
the  buckling  of  cylindrical  shells  subject  to  pure  bend- 
ing.  However,  these  classical  analyses  are  valid  only 
for  infinitesimal  deflection.   When  the  deflection  is  not 
very  small  and  has  a  magnitude  of  wall- thickness,  large 
deflection  theory  must  be  used,  and  the  high  order  terms 

of  deflection  are  to  be  considered.   In  general,  small 
deflection  theory  predicts  buckling  stresses  approximate- 
ly three  times  as  great  as  those  found  by  experiments. 
This  evidence  is  particularly  obvious  for  cylindrical 
shells  under  axial  compression.   The  discrepancy  between 
the  small  deflection  and  experimental  results  can  be 
explained  by  use  of  large  deformation  theory  which  was 
advanced  by  Donnell  (JLO)  and  von  Karman  and  Tsien  (jLl)  . 
In  Ref.  (j^) ,  Donnell  also  considered  the  initial  imper- 
fection in  the  shell  to  cause  lower  buckling  load.   A 
diamond- shaped  buckling  pattern  represented  by  a  function 


of  arbitrary  parameters  was  first  proposed  by  von  Karman 
and  Tsien  (11)  for  cylindrical  shell  under  axial 
compression.   Their  results  agreed  with  the  experimental 
results.   The  same  problem  was  also  solved  by  Kempner  (12) 
who  used  the  Ritz-method,  but  improved  the  von  Karman  and 
Tsien 's  solution  by  considering  five  free  parameters.   An 
investigation  of  the  buckling  of  cylindrical  shells  by  use 
of  Galerkin's  method  was  also  studied  by  Ekstrom  (13)  and 
Lu  and  Nash  (14) .   The  general  case  of  the  axially  compress- 
ed orthotropic  and  isotropic  cylinder  with  internal  pressure 
has  been  discussed  by  Thielemann  (15) .   The  general  equa- 
tions in  tensor  notation  for  isotropic  shells  were  derived 
by  Dill  (16)  who  also  included  the  effect  of  initial  imper- 
fections.  The  effect  of  initial  imperfection  on  the  buck- 
ling of  cylinders  under  different  loadings  was  studied 
in  Ref.  (_17 )  ,  (J^)  ,  (19),  (_20)  ,  (21)  and  (^2)  . 

The  general  instability  of  ring-stiffened  or 
stringer-stiffened  cylindrical  shell  by  the  use  of  small- 
deflection  relations  have  been  investigated  by  many  authors, 
Nash  (_23)  ,  Alfutov  (^4)  ,  McKenzie  (_25)  and  Kan  and 
Lipovskiy  (_26)  .   Further  study  in  general  instability  of 
ring-stiffened  cylindrical  shells  subject  to  external  hy- 
drostatic pressure  with  a  comparison  of  theory  and  experi- 
ment was  discussed  by  Galletly,  Slankard  and  Wenk  (22) • 


A  detailed  treatment  of  the  buckling  of  stiffened 
shells  has  been  discussed  by  Hedgepeth  and  Hall  (28) . 
The  problem  of  stability  of  orthotropic  shells  was  studied 
by  Becker  and  Gerard  (_29)  .   Further  study  in  general 
instability  of  orthogonally  stiffened  plates  or  cylindrical 
shells  has  been  discussed  by  Huffington,  Jr.  (30) ,  Neut 
(31)  and  Becker  (^2) . 

In  Ref.  (J2)  '  ^^   analysis  of  orthotropic  cylinders 
was  made  by  Taylor  whose  derivations  were  based  upon  the 
use  of  a  form  of  Donnell's  equation.   Recently,  an  analy- 
tical investigation  on  the  instability  modes  of  ortho- 
tropic cylinders  subject  to  compressive  or  bending  load  was 
made  by  Block  ( J3)  .   He  used  the  Dirac  delta  function  for 

rings  in  the  linear  equilibrium  equation  and  considered 
the  pressure  load. 

The  large  deflection  relations  applied  to  the  prob- 
lem of  the  instability  of  ring-stiffened  cylindrical  shells 
derived  energy  method  was  investigated  by  Lee  (3^) . 

Experimental  work  on  the  ring-stiffened  cylinders 
or  unstiffened  cylinders  under  different  loadings  has 

been  performed  by  MaCoy  (_3^)  ,  Peterson  and  Dow  (_36) , 
Lundquist  (J7)  ,  Harris  (_38)  ,  Fung  (29)  ,  and  Goree  (40)  . 


The  theory  of  anisotropic  elasticity  has  been  dis- 
cussed by  Hearmon  (41)  and  Lekhnitsky  (42) .   A  general 
discussion  on  the  theory  of  anisotropic  shells  has  been 
made  by  Ambartsumyan  (43) .   A  linear  solution  of  aniso- 
tropic cylinders  was  investigated  by  Cheng  and  Ho  (44) . 

The  purpose  of  this  study  is  to  investigate  the 
general  instability  of  cylindrical  shells  with  a  system 
of  stiffeners  which  is  mounted  at  various  angles.   The 
cylinder  is  internally  pressurized  and  under  one  of  the 
following  types  of  loading: 

a.  Uniform  axial   coirpression 

b.  Uniform  bending 

c.  Uniform  transverse  load. 

The  stiffeners  are  not  very  close  but  discretely  located 
and  thus  the  anisotropic  relation  does  not  apply.   The 
Kirchnoff  hypothesis  is  assumed;  the  eccentricity  of 
stiffeners  disregarded.   The  method  of  the  solution  is 
first  to  assume  the  deflection  function,  then  the  stress 
function  is  found  from  compatibility  equation  which  is 
carried  out  by  the  use  of  nonlinear  large  deflection  theory. 
The  expressions  for  the  energy  stored  in  the  shell  and 
stiffeners  as  well  as  the  work  due  to  the  external  loads 
are  obtained.   The  large  deflection  terms  and  the  effect 
of  initial  imperfections  in  the  strain  displacement 


equations  are  considered.   The  Ritz-method  is  used  to  find 
the  governing  equations  of  instability  of  stiffened  shells 

subject  to  these  three  types  of  loading  used.   Hence,  the 
minimum  loads  are  found  from  these  equations. 

Numerical  examples  are  worked  out  for  the  cylinders 
with  the  inclined  stiff eners  at  different  angles.   The  mini- 
mum strength  of  shell  versus  different  rigidities  of  stiffen- 
er  is  also  calculated.   The  relation  between  the  minimum 
strength  of  pressurized  shell  and  internal  pressure  is 
plotted  in  Fig.  15.   The  computation  was  carried  out  on  the 
IBM  709.   The  general  computer  programs  written  in 
Fortran  II  to  find  these  minimum  stresses  for  all  three 
types  of  loading  are  listed  in  Appendix  B.   Some  tests  have 
been  made  on  the  shell  with  inclined  stiffeners  at  angles 
of  Y  =  30°,  45°,  and  60°   under  axial  con^aression  and  bend- 
ing load.   Some  of  the  experimental  data  are  plotted  in 
Fig.  18  to  compare  with  the  theoretical  results. 


CHAPTER  II 

BASIC  RELATIONS 

2.1  Compatibility  Equations 

In  the  present  investigation,  the  following 
assumptions  are  made: 

a.  Elements  normal  to  the  unstrained  middle  sur- 
face of  the  shell  remain  normal  to  the  strained  middle 
surface. 

b.  The  material  follows  Hooke's  Law. 

Let  X  and  y  be  measured  in  the  axial  and  the  circum- 
ferential direction  in  the  median  surface  of  the  underformed 
cylindrical  shell  (Fig.  1),  w  equal  the  total  radial  deflec- 
tion, and  Wq   represent  the  initial  deflection  in  the 
radial  direction.   The  strain-displacement  relations  are 
(22): 


Tfy^    ,     X   C   ^ar\  ^    1  /  ^cj; 


^'  =  ^  +  i( 


2jx  '  2  V  ^x  /     Z     \>jc 


2 


""^    ^^    75X  -dX       ^^     ^X    -d^ 


8 


In  the  above  relations,  the  higher  order  terms  of  deriva- 
tives involving  u  and  v  are  neglected,  since  these  dis- 
placements are  small  compared  to  the  radial  displacement  w. 

The  stresses  and  the  strains  in  the  median  surface 
of  the  shell  in  the  case  of  plane  stress  are  related  to 
each  other  by  the  following  equations: 


^.= 


i'P' 


(ex-^-^6^) 


c  =       G 


(2) 


By  substituting  Eqs .  (1)  into  Eqs.  (2),  the  median  surface 
stress-displacement  relations  are  obtained: 


5x  -  1371(^7'^  2"  ^3x 


1    ,?cJ-.2    1   /^U7./  7jvll2f<^\ 


ZriT       l/>af 


_  1  (l}fl]'  _ 


,uf-~uT„\ 


1  f  ^f _  J.  fi^» 


&.  - 


■^xi, 


m+7)  L  ^y     ?x      "^x  ^^      ^x   21^  j 


(3) 


The  conditions  of  equilibrium  can  be  satisfied  by  using  the 
well-known  Airy  stress  function  F  which  is  defined  as 


(4) 


Eliminating  the  variables  u  and  v  in  Eqs.  (3)  and  (4),  the 
following  relation  between  stress  function  F  and  the  radial 
component  of  the  displacement  w  is  obtained. 


(5) 


This  equation  expresses  the  condition  of  compatibility  be- 
tween stress  and  strain.   It  was  first  obtained  in  the  pre- 
sent form  by  Donnell  (10)  . 

In  general,  Wq  is  unknown.   For  simplicity,  w  is 

assumed  to  be  proportional  to  w    Ref.  (10) ,  (17) ,  (20) , 
and  (22)  .   Thus  one  can  define 

■p  =  —    =   imperfection  ratio       (6) 

where  p  is  independent  of  x  and  y.   With  the  relation  from 
Eqs.  (5),  (6),  the  compatibility  equation  is  expressed  as: 


10 


(n^)-V  =  E(Hr).f(||-)- 


^ar    ^^or' 


MX^y/  15X}    Bryf 


(7) 


2.2   Energy  Expressions 

The  strain  energy  in  the  shell  and  stiffeners  and 
the  potentials  due  to  external  load  are  found  below. 

a.   The  extensional  strain  energy  in  the  shell  can 
be  written  as: 


IJe  = 


Z^E) 


:i 


rp   rv 


'dx^    ^y 


dx  dy 


(8) 


b.   The  bending  strain  energy  has  the  following 
expression: 

^  ?      .       «  \>Y^         >-r2>'  Wv2        ^-r2/    ^>h2        >h2/ 


2J 

+  ( 


Oj    0 


^X^       -bX^'    '^^^      arj2 


)^Zll-2)')[ 


-.1^(2^- 


Ir'oJ-^ 


a-xB5       3-x?^ 


n 


Ud. 


(9) 


where 


D  = 


Et 


c. 


The  potential  due  to  edge  bending  of  the  shell 


is 


11 


IL  = 


PL 
0  J 


ZTTR 


5. 


ixd^ 


(10) 


where  5^.  is  applied  peak  bending  stress  and 

lii__i  l^lL   _;pl!L^_  1  (— f+i  1— f 
Tix  "~  E  Um^    7^r^/   2  ^>r^  2  V75X/ 


?JX    E  \^^2 


d.   The  work  of  the  external  force  applied  at  the 
ends  of  the  shell  can  be  calculated  as  the  product  of  the 
applied  force  and  the  change  in  length  of  the  shell  as 
follow: 


Uc=6c* 


'i^ 


d.xd^ 


(11) 


where   f:  is  constant  stress  over  the  shell  thickness, 
^c 

e.   The  potential  due  to  the  internal  pressure  p  is 


U* 


ZTTR 


■f  [uT  —  ujh]    d.xd]j 


(12) 


The  change  in  volume  due  to  the  displacements  u,  v  are 
neglected  (11) . 

f.   The  potential  due  to  the  edge  bending  on  fixed 
end  by  transverse  shear  load  P  which  applied  at  the  free 
end  (Fig.  3)  is 
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g.   The  potential  due  to  transverse  shear  load  P 
applying  on  free  end  is 


IIs=- 


.  0  J  o 


7TR 


>y 


?>x 


)  dx 


+ 


jOi   0 


^uT   "^uT 


^x  :2.^ 


dxfl^ 


(14) 


where  P  is  average  transverse  shear  load. 

h.   The  bending  strain  energy  of  stiffeners.   The 
stif feners  are  assiimed  in  parallel  with  the  y  ,  y  coordi- 
nates lines  and  the  principal  directions  of  the  cylindrical 
shell  coincide  with  x,  y  lines  (Fig.  2) .   Let  the  subscript 
k  be  used  for  k^^  stiffener  which  is  inclined  at  an  angle 
Y-  with  the  generator  of  the  cylinders,  and  is  parallel 
with  y  -  line  and  normal  to  y'   -  line.   Thus  the  bending 
strain  energy  in  the  k   stiffener  is 


-h 


2.    ,  \2 
O 


(15) 
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where  N,  denotes  the  number  of  the  stiffeners  in  y. 
k  1 

-  direction,  and  E^  I^,  represent  Young's  modulus  and  the 
moment  of  inertia  of  the  k^^  stiffener,  respectively.   The 
eccentricity  due  to  stiffeners  is  neglected.   The  limit 

Ji.j^  is  the  length  of  the  stiffener  in  ^^  -  direction. 
Similarly,  the  bending  strain  energy  in  the  j^^  stiffener 
which  is  parallel  with  y,  -  line  and  normal  to  y'^^  -  line 

(Fig.  2)  is 


±L   C  T 
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(16) 


The  siobscript  j  is  used  for  j'^^  stiffener  which  is  inclined 

at  an  angle  of  iT  with  the  generator  of  the  cylinders. 

Where  N.  denotes  the  number  of  the  stiffeners  in  T ~ 
J  ^ 

-  direction,  E.,  I.,  represent  Young's  modulus  and  the 
moment  of  inertia  of  the  j^^  stiffener,  respectively.   The 
limit  I  .  is  the  length  of  the  stiffener  in  ^  j  "   direction. 

i.   The  torsion  strain  energy  of  the  k   and  j 
stiffeners  are 


^ 


Jo 


^^1 


y>o 


(17) 


14 
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'^^2~^^\       JV=0 


C^a; 


(18) 


where  G  and  J  represent  the  shear  modulus  and  the  polar 
moment  of  inertia  of  stiffeners,  respectively.   The  siob- 
scripts  j  is  for  stiffeners  in  "Jf   -  direction  and  k  in 
y   -  direction.   In  the  present  analysis,  the  inclined 
angles,  y   and  y  ,  are  considered  in  axial  symmetry 
(i.e.  TTj^  =  ^2  =  2'  ). 


CHAPTER  III 


INSTABILITY  OF  A  STIFFENED  CYLINDRICAL  SHELL  UNDER 
AXIAL  COMPRESSION  AND  INTERNAL  PRESSURE 


3.1   Deflection  Pattern  and  Stress  Function 

Let  a  deflection  shape  assume  the  form  (JL^) : 

W  =  b,+  b,COS^  COS^ik^COS'-^  +  h^COSiM   (19, 

where  m  and  n  are  the  numbers  of  waves  in  the  axial  and 
circumferential  directions,  respectively.   In  Eq.  (19), 
h-^   is  not  an  independent  parameter  but  is  used  to  satisfy 
the  condition  of  periodicity  of  circumferential  displace- 
ment (J^) .   It  is  found  later  to  be  a  function  of  bj,  b^ 
and  b^.   The  boundary  conditions  of  the  displacements  at 
the  ends  of  the  stiffened  cylinder  are  disregarded.   This 
simplification  is  justified  by  the  experimental  findings 
of  N.  Nojima  and  S.  Kanemitsu  as  discussed  in  Refs.  (11) 
and  (12) .   It  was  found  that  there  is  no  appreciable  length 
effect  when  the  length  of  the  cylindrical  shell  is  greater 
than  1.5  times  the  radius  of  the  shell.   Therefore,  the 
long  stiffened  cylindrical  shell  is  considered  hereon. 


15 
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The  corresponding  stress  function  is 

F  = -e^l+f  |-Va„cos^cos^+aa.cos^x 


The  coefficients  a,,,  ^22'  ^20'  ^02'  ^31  ^^'^   ^13 
are  found  in  term  of  b2/  b3  and  b^  from  the  compatibility 
equation  (22)  .   These  coefficients  are  found  to  be: 


(21) 


where 


(22) 
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Physical  considerations  lead  to  the  conclusion  that 
the  circumferential  displacement,  v,  must  be  a  periodic 
function  of  y.   This  is  expressed  in  the  condition  that 
^v/  ^y  contains  no  constant  terms  or  functions  of  x  alone, 
Using  the  second  of  Eq.  (1)  together  with  Hooke's  law,  the 
following  expression  for  ?v/  3y  can  be  obtained; 

By  substituting  Eqs.  (4),  (19),  and  (20)  into  Eq.  (23)  it 
is  found  that 


'd\r  ^  1  ;VK 


i^V6cj-(i-r^) 
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^^^\-V) 


.  (24) 

"T  terms  involving  periodic  functions 

Since  V  must  be  a  periodic  function  of  y,  the  constant 
term  on  the  right  side  of  Eq.  (24)  must  be  equal  to  zero. 
Thus, 

where 


t--^  K-Y^  (26) 
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3.2  Expressions  of  Total  Potential  Energy 

After  substituting  the  expressions  for  w  (19)  and 
F  (20)  into  the  relations  given  in  Chapter  II,  section 
2.2,  the  various  energy  terms  are  obtained  after  integra- 
tion.  They  are 
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(32) 
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h 


1^=-]^    '     C,  =C05^,     '     Ca^SiNr, 


p  =  -li- 
^j    L 


C3=C0S^2 


C4=SiNr, 


(35) 


The  expressions  R, ,  R2,  R3,  and  R-  are  very  lengthy  and  can 
be  neglected  if  the  stiffeners  are  located  uniformly  and 
not  far  apart.   They  are  listed  and  discussed  in  Appendix  A. 

The  total  potential  of  the  system  is  the  sum  of  the 
strain  energies  and  the  potential  of  the  applied  loads,  as 
follows: 


1^=-  ire+F,+Fctu;+ii/ii±T];.,-bir., 


(36) 


In  nondiraensional   form: 
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(37) 

3. 3  Minimization  of  Total  Potential  Energy 

The  variation  of  the  total  potential  with  respect 
to  each  of  the  arbitrary  parameters  must  vanish  for  equi- 
librium, this  yields 
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After  differentiation,  the  following  three  equations  are 
established: 
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For  brevity,  Eqs.  (39),  (40)  and  (41)  may  be 
rewritten  as: 
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Eliminating  b2  and  ^  froinEqs.  (42),  (43),  (44), 
the  following  equation  is  obtained: 
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Equation  (49)  is  the  governing  equation  for  the  ' 
stationary  load  of  a  stiffened  shell  under  axial  compression, 

The  numerical  calculations  have  been  carried  out  on 
the  IBM  709  and  the  computer  program  used  is  attached  as 
Appendix  B . 


CHAPTER  IV 


INSTABILITY  OF  A  STIFFENED  CYLINDRICAL  SHELL 
UNDER  BENDING  AND  INTERNAL  PRESSURE 


4.1  Deflection  Pattern  and  Approximate  Stress  Function 

When  a  cylindrical  shell  is  subjected  to  either 
pure  bending  or  eccentrically  applied  compression,  an 
approximate  form  of  the  deflection  pattern  is  assumed 

(14): 

^  =  b,  +  C05Uj(b,COS^C05^+b3C0S^^+b,C0S^) 


(53) 


This  approximate  deflection  pattern  is  obtained  by 

multiplying  the  pattern  of  the  axial  compression  by 

2 
cos   (y/2R)  which  signifies  the  localized  buckling  on  one 

side  of  the  stiffened  shells.   The  effects  of  the  end 

restraints  are  also  neglected. 

The  corresponding  stress  function  F  is  proposed 
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+  a„  COS  ^  COS  f  +  a„  cos  f  cos^ 


(54) 


The  stresses  ^  anci  g"   are  the  average  axial 
compression  and  peak  bending  stress,  respectively,  and  are 
positive  for  compression. 

The  coefficients  a^o'  ^02'  ^11'  ^22'  ^31  ^^^   ^13  ^^ 
Eq.  (54)  can  be  expressed  in  terms  of  b  ,  b   and  b  by  the 
Galerkin  method.   This  method  establishes  the  following 
set  of  equations: 
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where 
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After  the  expressions  for  uT  in  Eq.  (53)  and  F  in  Eq.  (54) 
are  substituted  in  Eq.  (55)  and  the  integration  carried 
out,  the  following  relations  are  found: 
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where  p  -  yU  .  ^  and  b^^  are  same  as  given  in  Eqs.  (6)  and 
(22),  respectively. 

In  order  to  find  the  b,  as  a  function  of  b2,  h^   and 
b^,  we  derive  as  in  section  3.1  the  relation: 
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(58) 

4.2   Expressions  of  Total  Potential  Energy 

The  total  potential  of  the  system  is  the  s\im  of  the 
strain  energies  and  the  potential  of  the  applied  loads.   By 
adding  the  Eqs.  (8),  (9),  (10),  (11),  (12),  (15),  (16),  (17) 
and  (18),  the  total  potential  may  be  expressed  as  follows: 

(59) 

If  w,  F  and  b,  as  given  by  Eqs.  (53),  (54)  and  (58) 
respectively  are  sxibstituted  in  Eq.  (59)  and  the  integration 
carried  out,  the  total  potential  is  obtained  in  the  follow- 
ing form: 
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where  "r  ,  R  ,  R  ,  R  are  also  the  portion  of  the  energy  of 
the  stiff eners  which  have  the  similar  expressions  as  R^^, 
R,,  R3  and  R. .   Since  they  are  neglected  in  the  numerical 
calculation,  their  detailed  expressions  will  not  be 
elaborated  here.   From  the  discussion  in  Ref.  (14)  it  is 
indicated  that  m  has  a  greater  magnitude  when  R/t  is 
greater.   In  the  present  analysis,  we  are  concerned  only 
with  extremely  thin  shells?  hence  this  ratio  is  large. 
Therefore,  for  practical  purposes  l/m^  and  1/m'*  are 
negligible  compared  to  unity.   Thus,  from  Eq.  (60),  the 
total  potential  can  be  simplified  as  follows: 
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where 
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4.3  Minimization  of  Total  Potential  Energy 

The  total  potential  energy  must  be  a  minimum  when 
the  structure  is  in  equilibrium  and  this  condition  leads  to 
the  following  equations 


(63) 


Substitution  of  Eq.  (61)  into  Eq.  (63)  leads  to 
three  equilibrium  conditions  as  follows: 
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Eliminating  R    and  "Bj  from  Eqs .  (64),  (65)  and 
(66),  the  following  equation  is  obtained: 
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Eq.  (70)  is  the  equation  for  the  stationary  load  of 
a  stiffened  shell  under  bending  load.   Numerical  calcula- 
tions of  the  minimum  stress  parameter  $2  ^^^  given  in 
Chapter  VI. 


CHAPTER  V 


INSTABILITY  OF  A  STIFFENED  CYLINDRICAL  SHELL  UNDER 
TRANSVERSE  LOAD  AND  INTERNAL  PRESSURE 


5.1   Deflection  Pattern  and  Approximate  Stress  Function 

When  a  long  cylindrical  shell  is  sxibjected  to  trans- 
verse shear  load  at  one  end,  an  approximate  form  of  the 
deflection  pattern  can  be  assumed  as  in  the  bending  case 
in  the  previous  chapter  in  Eq.  (53).   The  deflection 
shape  is 

Corresponding  to  this  deflection  pattern,  the  stress 
function  F  is  proposed 

+  a,,  cos^f^  cos  ^  +  a,„  cos^  +  a,,w5^ 
+  a3,c..^cos2^«+a„cos2^cosi^ 

(74) 
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The  relations  between  the  coefficients  aTi.  ^22'  ^02' 
^20'  ^13'  ^31'  ^^^   ^2'  ^3'  ^4  ^^®  given  in  Eq,  (57), 
and  15^   is  the  same  as  given  in  Eq.  (58). 
5.2  Expressions  of  Total  Potential  Energy 

The  total  potential  energy  of  this  system  is  the 
sum  of  the  strain  energy  and  the  potential  of  the  applied 
loads. 


TL==lf.^V:+V'^^TJs^JL^'UU^.,t^r,^^r.-,      (7 


(75) 


If  w,  F,  and  bj^  as  given  by  Eqs.  (53),  (58),  and 
(74),  respectively,  are  substituted  in  Eq.  (75),  and  the 
integration  carried  out,  the  total  potential  is  obtained 
in  the  following  form: 


11  = 


UbR 


T?^2  —2 


=  -p  +  V-  2(H-j;3P>  (l-h^;(fj  p 


(l-r)(¥^^^)^  +  f-2^l^ 


-2  (i+M.^f 


I 2 


-t 


-H  «/ 


4  ^^ 
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(i-vf  r-T- 


-h 


t: 


2  L 


2  ,   6^^  ,2,1 


I'^^^^d-^i^^j 


U(i4^'-J^ 


(i+^'r  + 


4^,  4- 


(76) 


where       ,  ,   IT.   I  I    I  I  .     are  the  same  as  those 
obtained  in  the  bending  case  Eq.  (60) . 


wh 
and 


ere  fl,  ,     n.   and  y    are  given  in  Eqs.  (22)  and  (6), 

PL 


P 


Ei\Hr 


{ID 


For  the  same  reason  as  mentioned  in  the  previous 
chapter,  l/m^  and  1/m^  can  be  neglected  compared  to  unity 
in  Eq.  (76) .   The  total  potential  energy  is  finally  ob- 
tained in  the  following  form: 


U, 


5  ^i 


Ka2  7^2 


P+t'-2(l  +  ^)P'+(l-Hi^J(f3  P 


—  2 


-(l-rK|^+b^'}:^+2t 


4f 


—  2 


- 1 


4  ^" 


4^' 
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+ 


ii-ry 


192  ^^{I-J)^J 


—  2 


(78) 


where   a^ii  /    ^22'    '^02'    ^20'    ^31'    ^^^  ^13   ^^®  given   in 

Eq.     (62). 

5.3  Minimization  of  Total  Potential  Energy 

By  variation  of  the  total  potential  energy  in 
Eq.  (78)  with  respect  to  each  of  the  arbitrary  parameters, 
three  simultaneous  equations  are  obtained.   From  these 
three  equations,  two  arbitrary  parameters  (b2  and  A  )    can 
be  eliminated.   Finally,  a  governing  equation  f or . the 
stationary  load  of  a  stiffened  shell  under  transverse  shear 
load  is  obtained: 


M,5'h-M2  5,+  M3  =0 


(79) 


where 


3  Z 


(80) 


M-,  ,    M2,    and  Mo    are   given   in  Eq.     (71) 


CHAPTER  VI 


NUMERICAL  EXAMPLES 


The  minimum  stresses  are  to  be  found  from  Eqs.  (49), 
(70),  or  (79)  for  the  various  types  of  loading.   In  each 
case,  the  stress  parameters  (  ^    ,    ^^  ^^^   X3^  ^^®  func- 
tions of  three  free  parameters  07  ,  "^i*  ^^^  y^*      While  xt 
is  the  wave-length  ratio,  the  expression  'Y!    and  %,,    can 
be  considered  as  deflection  parameters  ratio.   The  minimum 
values  are  found  numerically  at  various  inclined  angles, 
stiff eners  rigidity  and  the  internal  pressure.   Numerical 
calculations  were  made  on  the  IBM  709.   Program  I  is  for 
finding  minimum  stress  ^   for  axial  compression  case. 
Program  II  is  for  finding  minimum  stress  parameters,  ^2 
and  5 3.   The  minimum  stress  can  be  found  with  any  given 
numerical  data.   In  the  following,  some  examples  are  given. 

6.1  Cylinder  under  Bending  and  Internal  Pressure  or 
Transverse  Load  and  Internal  Pressure 

In  this  part  of  calculation,  the  cylinder  is  assumed 

perfect,  i.e.,  p  =  0  [Eq.  (6)]  ,  and  the  angles  of  the 

inclined  stiffeners  with  the  generator  are  equal,  i.e.. 
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-^      =    y„  =  y  .   Poisson'B  ratio  is  1/3.   The  general 
characteristics  are  assumed: 


Mi  Na 

■>c —  -t     —  —  ^  "  I  —  — 


(81) 


2 


By  substituting  Eq.  (81)  into  Eq.  (70),  the  value  of  ^ 
as  a  function  of  Tf  ,    'yt      and  yW.  is  found.   First,  ^ 
versus  77  is  plotted  for  various  values  of  yU  at  a  fixed 
value  of  7^  ^,  in  Figs.  4,  5,  6,  7,  8,  and  9.   Minimum  ^2 
found  from  each  of  these  curves  is  called  $2'  -?/  •  ^^ 
should  be  equivalent  to  the  value  found  from  the  relation 
^$2/^^  ~  ^'      Ths"^  ?2  y1    versus^  is  plotted  for 

various  71       in  Fig.  10.   Minimum  J  9  v  found  from  each  of 
-'■  ^1  % 

these  curves  is  called  5  2  77  >u   Plotted  in  Fig.  11.   It 
should  be  equivalent  to  the  value  found  from  the  relation 
^  ?2  7/  /^'•^  ~  ^*   This  minimum  value  is  the  dimensionless 
critical  stress  (  ^     +  3/2  5"  )    for  the  shell  with  in*?* 
clined  stiffeners  at  an  angle  of  ^  =  45*^  with  the 
generator.   By  the  same  procedure  as  in  the  previous  case. 
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the  curves  of  $  ^  ^     versus  71   ,  for  various  inclined 
^  2,yi  ,M-  ^  1 

stiffeners  at  an  angle  of   JT  =  0°  (stringers),  30°,  60°, 
and  90°  (rings)  are  found  and  shown  in  Fig.  12.   The 
minimvun  value  of  each  of  these  curves  is  plotted  in 
Fig.  13.   The  relation  between  rainimum  stress  parameter 
S^  ^     and  stiffeners  rigidity,  E  I  Ji  ,  is  shown  in 
Fig.  14.   The  relation  between  minimum  stress  parameter 

1E>  ^    „     and  internal  pressure  p  is  also  plotted  in 
-J-  2,  7^  ,  M- 

Fig.  15. 

If  ^   is  replaced  by  ^   in  Eq.  (79),  the  equation 
is  the  same  as  the  one  obtained  in  the  bending  case. 
Therefore,  the  minimum  stress  for  the  case  of  transverse 
shear  can  be  found  by  using  '^  instead  of  5  o  >f 

in  all  the  figures  from  Fig.  4  to  Fig.  15.   Where 


=  1-^1 


cr. 


6.2  Cylinder  under  Axial  Compression  and  Internal 
Pressure 

By  the  same  procedure  as  in  the  previous  section  we 

can  find  the  minimum  stress,  fT   from  Eq.  (49)  by  Program 

c 

I  in  Appendix  B.   In  this  part  of  calculation,  a  numerical 
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example  is  given  for  comparison  of  the  effect  of  the 
Imperfection  In  shells  at  various  Inclined  angles. 
This  is  plotted  in  Fig.  16. 


CHAPTER  VII 


EXPERIMENTAL  INVESTIGATION 


A  group  of  tests  for  the  general  instability  of  v 
rings  and/or  stringers  stiffened  shells  have  been  made  by 
many  authors,  but  up  to  the  present  time,  no  work  has  been 
done  on  the  experimental  investigation  of  the  general 
instability  of  inclined  stiffened  cylinders  subject  to 
axial  compression  or  bending.   Therefore,  it  was  necessary 
to  conduct  a  series  of  tests  in  order  to  compare  experi- 
mental with  the  theoretical  results. 
7.1  Models 

The  models  used  in  these  series  of  tests  were 
stiffened  cylindrical  shells  constructed  from  Du  Pont 
Mylar  of  1000  gage  (0.01  inches),  type  A.   Tests  indicated 
that  the  Mylar  sheet  has  a  Young ' s  modulus  E  varied  from 
550,000  psi  to  780,000  psi  and  Poisson's  ratio  of  1/3. 
In  the  numerical  calculation  to  follow,  the  value  of  E 
is  taken  to  be  700,000  psi. 

All  of  the  models  were  made  by  rolling  the  Mylar 
sheet  around  a  thick-walled  steel  tube  and  Joining  the 
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ends  with  a  3/4  inch  wide  strip  of  double-faced  Scotch 
tape.   It  appears  that  these  joints  might  stiffen  that 
part  of  the  cylinder  appreciably.   However,  reports  from 
various  investigators  indicate  that  this  is  not  the  case, 
since  buckling  waves  appear  across  the  joint  with  no 
noticeable  change  in  pattern.   The  present  tests  have 
confirmed  this  observation.   The  inside  radius  of  the 
cylinders  was  4  inches,  while  the  length  was  9  inches. 
The  stiffeners  were  made  by  cutting  the  Mylar  tape  to  a 
band-width  of  3/8  inch.   The  thickness  of  each  layer  of 
the  tape  was  0.01  inch.   All  the  specimens  had  stiffeners 
made  of  two  layers  of  the  Mylar  tape  (i.e.,  nominal  thick- 
ness of  the  stiffeners  =  0.02  inch).   The  double-faced 
Scotch  tape  was  used  to  join  the  layers  of  f^lar  tape.   It 
was  also  used  to  bond  the  stiffeners  with  the  cylindrical 
shells.   The  adhensive  tape  was  not  effective  when  more 
than  two  layers  of  Mylar  tape  were  used  as  stiffeners. 
The  spacing  between  two  neighboring  stiffeners  was 
1  1/4  inches.   The  stiffeners  were  inclined  at  angles 
of  'J^'  =  30°,  45°,  and  60°  with  the  generator  of  the 
cylinders.   In  order  to  provide  additional  information 
for  the  numerical  example  of  the  previous  chapter,  the 
dimensions  of  the  stiffened  shells  used  in  the  experiments 
were  the  same  as  those  considered  in  that  exait^jle. 
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7.2   Test  Result 

The  cylinders  were  mounted  vertically  on  the  test 
machine.   The  upper  adaptor  has  two  pivot  pins  fixed 
diametrically  at  its  edge.   The  pins  are  supported  by 
horizontal  bearings  which  are  fixed  to  the  frame  of  the 
test  stand.   When  a  moment  is  applied  to  make  the  upper 
adaptor  rotate  about  its  pivot  pins,  this  transmits  the 
bending  moment  to  the  cylinder.   The  lower  adaptor  is 
connected  to  a  circular  plate  which  can  be  moved  up  and 
down  so  that  the  axial  loading  can  be  transmitted  to  the 
cylinder.   Bending  load  is  applied  by  means  of  a  lead 
screw  which  pulls  a  cable  up  through  a  bearing  which  is 
fixed  on  the  frame  of  the  machine. 

The  bending  moment  and  the  axial  compressive  loading 
were  found  after  calibrating  the  readings  from  a  strain 
gage  indicator.   The  test  results  of  the  stiffened  shells 
subject  to  axial  compression  or  bending  without  internal 
pressure  at  an  angle  of  if  =  30°,  45°,  or  60°  are  tabu- 
lated in  Table  I.   Some  of  the  results  are  plotted  in 
Fig.  18.   The  deformed  pattern  after  buckling  at  p  =  0 
are  shown  in  Fig.  17.   The  buckle  pattern  formed  across 
the  seam  with  no  distortion  in  shape  and  none  of  the  seams 
failed  during  buckling. 


CHAPTER  VIII 
DISCUSSIONS  AND  CONCLUSIONS 

In  the  previous  chapters,  nonlinear  analysis  on 
cylinders  with  inclined  stiffeners  have  been  made  by  the 
energy  method.   From  the  numerical  examples  the  following 
observations  are  made: 

1.   From  Fig.  13,  it  can  be  observed  that  the  mini- 
mum stress  parameter  5  o  ^  j,     increases  with  increasing 
stiffeners'  rigidity,  and  the  optimum  inclined  angle 
varies  with  the  rigidity.   In  the  present  example,  it  has 
been  found  that  the  most  effective  inclined  angle  is  in  the 
neighborhood  of  ^  =  60°. 

2.   In  Fig.  14,  strength  of  cylinder  with  stiff ener 
inclined  at  /  =  45°  is  compared  with  that  of  ring- 
stiffened  cylinder.   At  smaller  stiffeners'  rigidity,  the 
ring-stiffened  cylinder  is  stronger,  but  at  higher  rigidi- 
ty, the  inclined  stiffened  cylinders  have  more  strength. 
The  minimum  stress  of  ring-stiffened  shells  approaches 
to  the  buckling  stress  found  from  small-deflection 
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solutions  as  a  limit,  while  the  inclined  stiffened 
cylinder  continues  to  increase  with  increasing  rigidity 
(Eli).   In  the  same  figure,  the  consideration  of 
torsional  rigidity  (for  example  K  =  1/2)  is  compared  with 
the  result  including  the  bending  rigidity  only  (i.e., 
K  =  0) .   The  difference  between  these  results  increases 
with  rigidity,  E  I  J  .   However,  the  general  relation  of  the 
minimum  stress  with  other  parameters  will  remain  the  same. 

3.  Example  of  the  variation  of  $^  ^     versus  p 
is  shown  in  Fig.  15.   For  other  combinations  of  "^T,  K  and 

E  I  fi ,  the  curve  will  be  similar  but  different  numerically. 

4.  The  minimum  stress  of  cylinder  under  cixial 
compression  versus  the  various  inclined  angles  is  shown  in 
Fig.  16.   Imperfect  cylinder  (  P  =  0.3)  has  lower  strength 
as  expected.   The  evaluating  of  P  can  be  referred  to 
Ref.  (22).   The  effect  of  imperfection  for  stiffened' 
cylinders  is  smaller  in  comparison  with  the  unstiffened 
cylinder. 

5.  The  buckling  patterns  of  the  unpressurized- 
stiffened  shells  subject  to  axial  compression  are  shown 
in  Fig.  17.   The  deformed  pattern  after  buckling  is 
diamond- shaped  and  across  the  inclined  stiffeners  as 
expected. 
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6.   In  Fig.  18,  it  can  be  observed  that  the  theory 
for  the  bending  case  is  in  reasonable  agreement  with  tests 
on  the  Mylar  inclined-stiffened  cylindrical  shells. 

This  study  has  presented  the  approximate  solutions 
by  energy  method.   Nevertheless,  the  results  of  this 
analysis  should  give  some  insight  into  the  problems  of 
the  general  instability  of  stiffened  shells  with  inclined 
stif feners. 
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TABLE  I 


TESTS  OF  STIFFENED  THIN  MYLAR  CYLINDERS  UNDER  AXIAL 
COMPRESSION  OR  BENDING  WITHOUT  INTERNAL  PRESSURE 


Cylinder  Radius  =  4  inches 
Mylar  Nominal  Thickness  of  Shell, 


t  =   0.0075  inches 


Thickness  of  Stiffeners, 


=   0.02  inches 


Inclined 

Dimensionless 

Dimensionless 

Dimensionless 

angles 

rigidity  of 

cixial 

bending 

stif feners 

compression 

stress 

E  I  i 

?c 

5b 

45° 

0.059 

0.267 

0.306 

0.274 

0.317 

0.117 

0.282 

0.312 

0.294 

0.329 

0.308 

0.317 

0.302 

0.306 

0.47 

0.287 

0.321 

0.339 

0.382 

0.322 

0.323 

0.94 

0.368 

0.435 

0.328 

0.388 

0.335 

0.412 

0.336 

0.376 

1.47 

0.333 

0.365 

0.335 

0.365 

0.339 

0.358 

30° 

0.124 

0.281 

0.306 

0.315 

0.321 

0.992 

0.302 

0.353 

0.322 

0.335 

0.314 

0.353 

60° 

0.134 

0.288 

0.388 

0.322 

0.322 

1.072 

0.302 

0.353 

0.322 

0.370 

0.322 

0.370 
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X,  u 


Pig.  1.   Coordinates  and  Displacement  Components  of 
a  Point  on  the  Middle-Surface  of  the  Shell, 
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//  ^  \ 


Fig.  2.   The  Coordinate  System  of  the  Stiffened 
Shells  and  Stiffeners. 


(a)   Under  Axial  Compression  P 
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M 


M 


(b)   Under  Pure  Bending  M 


(c)   Under  Transverse  Shear  Load  P 


Fig.  3.   The  Dimensions  of  the  Stiffened  Shells  under 
Different  Loadings. 
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Comparison  of  Effect  of  the  Imperfection  Ratio 
on  the  Stiffened  Shells  Under  Axial  Compression 
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APPENDIX  A 


APPENDIX  A 


ENERGY  R^  IN  STIFFENERS 


To  neglect  the  energy  terms  R^ ,  R^,  R.  and  R^  in 
Eq.  (37),  we  may  refer  to  Ref.  (34).   In  that  study,  the 
numerical  results  have  indicated  that  the  minimum  stress 
is  only  slightly  affected  by  these  terms.   Also,  it  is  to 
be  noted  that  the  minimum  stress  is  nearly  independent  of 
the  locations  of  the  stiffeners.   In  the  present  analysis 
and  experimental  results,  it  indicates  that  the  shells 
buckle  in  multiple  waves  pattern.   This  implies  that  m^ 
and  n2  are  much  greater  than  unity.   Moreover,  from  the 
experiment  the  buckling  stress  is  not  affected  by  small 
change  in  the  stiffeners'  locations.   Thus,  the  summations 
of  sine  functions  of  the  energy  of  the  stiffeners  may  be 
assumed  very  small  when  compared  with  the  constant  terms. 
The  approximate  solution  can  be  much  sin^lified  by 
neglecting  these  R^  terms  (i  =  1,  2,  3,  4).   However,  the 
expressions  for  R^^,  R^,  R^  and  R  are  expressed  below  for 
the  purpose  of  reference. 
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L(7nCi.-37?C2} 


^£/z  (;/Z  C,-  3;?<2  .)PS  -•■  ,— 


^    «        U^r2C^^nC-^) 


9      r~ 


7"     T      ..-2  _  2 


+  Uii^^;^)  .«(«c,+.c,)(A|,,6t3b,«'c/c; 


R 


X  sm(e  7.c^-i^C2)(-^H,,^,,%e,;^-^-^  ^(^^.-^^.)(^;] 


4  li[(C,^^C,j"  ^,(,,^^^,^^    5^^   z(^c,..C,)^(Aj^ 
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^^r^^^^'  ZL(Z,-nc,)  '^'  ^(-V^^P(#}>2b3u\    ^^/ 


/c*  \   1    ■?  u     n      — — 


I 


xs<«^ncj(^)  +  2b,C,  -j^  stn(4,ic,i^}^. 


r  — 2 


Rr"^^  f^  lili  ^  (C3  +  ^  C,)^fc3-^  c,)  (^  S<.  Z^C,^} 


■'      ^-     pt     ^      -4 


+ 


R 


^LnC^ 


5m  2nc^  (^J  +  2 b2  1^3  u'iC.^uC^f  c/  f^ 


7^^3-72  0^) 


^^^^(-^^3-^C,)(^H,^^,^%,,^^    3.<-C3+3^C^)(Aj-| 
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5Ln{3>mc^^nc^){-^) 
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L(3WC3--^C4) 
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li3mC)~nc^) 


3a^  (3?nC^-nc^)(^)  -f- 


LirrfC^-\.nc^) 


-  \(:^M"b^b^C^ci 


^Liwc^+nc^) 
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y  5^  {'mc^-hnc4.)(^)j 


1^  ZL(7nC^-7iC^)  ^^K^j 


APPENDIX   B 


APPENDIX  8 
COMPUTER  PROGRAM  I 
STIFFENED  CYLINDRICAL  SHELLS  UNDER 

AXIAL  COMPRESSION  I 

1  READ  INPUT  TAPE  5,  10 , E , G, C, D, Q , R, H, P, I  I , I D, I F , J  I , JD , J  ' 

1F,KI,KD,KF 
10  FORMAT  (8F5.3,9I4) 

DO  50  K=KI,KF,KD  I 

Z=.Ol»FLCATF(K)  i 

DO  50  I-II,IF,ID 

X=.01»FLCATF{I  ) 

DO  50  J=JI,JF,JD 

Y=.G01*FL0ArFlJ) 

YS=Y*»2 

XS=X»»2 

XQ=XS»»2 

CS=C**2 

CQ=CS»»2  I 

US=D»*2  I 

DQ=DS»*2 

QS=Q*»2  I 

GQ=QS»»2 

RS=R»»2  \ 

RQ=RS»»2 

HM=1.-H 

HP=1.+H  i 

HMS=1./HM»*2 

S1=E»{CG+GG+XQ»(DQ+RQ)+6.»XS»(CS»DS+QS*RS)+G»( ( (C+X»D) 
1*(D-X*C) )*»2  +  ( {C-X»D)»(D*X»C) )»*2+( { Q  +  X»R ) • { R-X»Q  )  ) •»2 
2+( (Q-X*R)*(R+X*e) )»*2)/2.)/4.  i 

52  =  E»XQ»(Dt;  +  RG  +  G*(CS»DS  +  CS»RS)  )  | 

53  =  E*(CG+GG  +  G»(CS»DS+QS«RS)  )  j 
G1=1./(1.+XS)«*2 

G2=l./(9.+XS)**2 

G3=l./( l.+9.»XS)»*2 

AI  =  3.»{  l.+XS)»»2/32.  +  Sl 

A3=-{2.»{2.+H)»( l.+Z)»Gl+Z/2.)»XS 

A4=HP*A.*XQ*(G1»(1.+Z)»»2+G2»Z**2+G3)  \ 

A5=HP»( l.+XQ)/16.  I 

Bl=3.»XG/8.+S2 

B4=HP»2. »XQ»G1 »Z»»2 

B5=HP»XG*.5»(G3+(1.+Z)*G1) 

B6=-XS».25»G1 

CI=3./8.+S3 

C4=HP»2.»XQ*G1  j 

C5=HP»XQ».5*{G2+(1.+1./Z)*G1)  ' 

C6=-XS»(Gl+(l.+H)/8.)».25/Z 

Hl=C5+C6/Y-A5 

H2=B5+B6/Y 

D2=Al*H2-Bl*A5 

D3=(G1+A3»Y+A4»YS)»H2-B4»A5»YS 

F2=Al»(Cb+C6/Y)-Cl«A5  | 

F3=(G1+A3»Y+A4»YS)»(C5+C6/Y)-A5»(.25+C4»YS) 
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60 


70 


80 
20 
50 


1F3)» 

W2  =  - 
i3-(2 
2D2)/ 
hi3=H 
1»2/F 
23)»» 
DISC 
IF  { 
XIR= 
X2R= 
GO  T 
S=SQ 
X1R  = 
X2R  = 
GO  T 
WRIT 
FORM 
CONT 
GO  T 
END 


D2»D 

HMS 

HMS» 

.»H2 

F3) 

MS*2 

3-(H 

2-2. 

=  W2» 

DISC 

-W2/ 

XIR 

C  80 

RTF( 

(-ifJ2 

(-W2 

0  80 

E  OU 

AT  ( 

INUE 

0  I 


3»Hl*HI/F3»»2+F2*H2»H2/F3-H2*HI»(F2»D3/F3+D2)/ 

3»XS«P»(D2»D3»HI»(Hl+A5)/F3»»2+F2*H2»(H2-A5)/F 
»{H1+A5)*A5»A5-H2»A5-(H1+A5)»A5)».5»{F2»D3/F3+ 

.25»XQ»P«P»(D2«D3*(H1+A5)»»2/F3**2+F2»(H2-A5)» 

1  +  A'>)»(H2-A5)«{F2*D3/F3+D3)/F3)+D2»D2+(D3»F2/F 

♦D2*D3»F2/F3 

»2-4.«WI»W3 

)  30,  60,  70 

(2. 'WD 


DISC) 

+S)/(2.*W1) 

-S)/(2.*W1) 


TPUr  TAPE  6, 
10F8.3) 


20,  Z,X,Y,XlR,X2R,P,Wl,W2,W3fE 
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Symbols  Used  In  Computer  Program  I 

E  =  E  "l  X 
C   =  C^^ 
Q  =  C3 

H   =   r 

Y  =  7^ 


G 

= 

K 

D 

= 

S 

R 

= 

^=4 

P 

= 

P 

X 

= 

^ 

X2R 

= 

Cc 
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COMPUTER  PROGRAM  II 
SriFFENEO  CYLINDRICAL  SHELLS 
UNDER  BENDING  OR  TRANSVERSE  SHEAR 
I  READ  INPUT  TAPE  5,  10, £ , G, C , D, Q , R, H, P, I  I , I D, I F , J  I , JD , J 
lF,Ki,KD,KF  ! 

10  FORMAT  (8Fb. 3,914) 

00  50  K=KI,KF,KD  j 

Z=.01*FL0ATF(K) 

DO  50  I=II, IF, ID  ! 

X=.01»FLUATF(I)  i 

DO  50  J=JI,JF,JD  J 

Y=.001«FLUAIF( J) 

ZS=Z**2  1 

yS=Y»«2 

XS=X»*2 

XQ=XS»*2  I 

CS=C»«2  i 

CQ=CS»*2  j 

DS=D»»2  I 

DQ=DS»*2  I 

QS=U»»2 

QQ=US»*2  j 

RS=R»»2 

RQ=RS»»2  ; 

HM=l.-H  j 

HP=1.+H  ! 

HMS=1./HM*»2  1 

Sl  =  3.»E»(Cg*Q(J  +  Xg»{DQ  +  RQ)+6.»XS*(CS*DS  +  QS»RS)+G»{  ((C  +  X 
1»D)»(D-X*C) )»*2+( (C-X»D)*(D+X»C) )»»2+( ( U+X*R ) ♦ I R-X»Q ) ) 
^»*2*-{  (Q-X*R)»{R  +  X»Q)  )*»2)/2.J/8.  | 

S2=1.5*E»Xg»(DQ+RQ+b*(CS»DS+QS»RS) ) 

S3=1.5»E»(CQ+Qg+G*(CS»DStQS»RS) )  j 

Gl=l./( l.+XS)»*2 

G2=l./(9.+XS)»»2 

G3=l./(l.+y.»XS)*»2 

Al=9.»l l.+XS)»»2/64.+Sl  ' 

A3=-(1.5»(2.+H)«( l.*Z)»Gl+3.*Z/8.)»XS 

A4=HP»9.»Xy»(Gl»( l.+Z)»»2+G2»ZS+G3)/4. 

A5=HP»9.«( l.+XQ)/256.  ) 

Bl=9.»Xg/16.+S2 

B4=HP»9.»XG»Gl»ZS/3. 

B5=HP»9.»Xg»(( l.+Z)»Gl+G3)/32. 

B6=-(3.»XS*Gl)/i6. 

CI=9./16.+S3 

C4=HP»9.«XQ»Gl/8.  , 

C5=HP»9.*XU*( ( l.*l./Z)»Gl+G2)/32.  | 

C6=-3.»XS»(HP/8.+Gl)/(  16.»Z)  '] 

H1=C5+C6/Y-A3  I 

H2=B5+B6/Y  ; 

D2=A1»H2-B1»A5 

03=IG1+A3*Y+A4»YS)»H2-B4*A5»YS 

F2=A1»(C5+C6/Y)-CI»A5 
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F3=(G1+A3*Y+A4»YS)»{C5+C6/Y)-A5»(.25+C4»YS) 

W1=HMS*.25»{ D2»D3»HI»H1/F3»»2+F2»H2»H2/F3-H1*H2»(F2»D3 
l/F3+02)/F3) 

W2=-HMS*1.5*XS»P»(D2*D3»HI»(H1+A5)/F3»*2+F2»H2»(H2-A5) 
i/F3-(2.»H2»(Hl+A5)+A5»A5-H2»A5-(Hl+A5)*A5)».5»(F2*D3/F 
23+D2)/F3) 

W3=HMS*^.25*XW*P»P»{U2*D3»(H1+A5)»»2/F3*»2+F2*{H2-A5)» 
i»2/F3-(Hi  +  Al3)»(H2-Ab)*tF2»D3/F3  +  D3)/F3)+D2*D2+(D3*F2/F 
23)»*2-2.*D2*D3»F2/F3 

DISC=W2**2-4.»W1»W3 

IF  (DISC)  bOi    60,  70 
60  X1R=-W2/(2.»»J1) 

X2R=XiR 

GO  TO  80 
70  S=SQRTF(UISC) 

X1R=1-W2+S)/{2.»W1) 

X2R=(-W2-S)/(2.»W1) 

GO  ro  80 

60  WRITE  OUTPUT  TAPE  6, 
20  FORMAT  (10F8.3) 
50    CONTINUE 

GO  TO  1 

END 


20,  Z,X,Y,X1R,X2R,P,W1,W2,W3,E 
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Symbols  Used  In  Computer  Program  II 


E 

=      E    I 

I 

G 

=      K 

C 

=      ^1 

D 

=      S 

Q 

=      C3 

R 

=      C4 

H 

=  r 

P 

=      P 

Z 

=  ^, 

X 

=  >t 

Y 

=  ^ 

X2R 

=      2 

_  =  2(  5"^  +  3/2  &)  or 

(X2R  =   2  ^    =  "p   for  transverse 
load) 
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